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We present atom-interferometer tests of the local Lorentz invariance of post-Newtonian gravity. 
An experiment probing for anomalous vertical gravity on Earth, which has already been performed 
by us, uses the highest-resolution atomic gravimeter so far. The influence of Lorentz violation 
in electrodynamics is also taken into account, resulting in combined bounds on Lorentz violation 
in gravity and electrodynamics. Expressed within the standard model extension or Nordtvedt's 
anisotropic universe model, we limit twelve linear combinations of seven coefficients for Lorentz 
violation at the part per billion level, from which we derive hmits on six coefficients (and seven 
when taking into account additional data from lunar laser ranging). We also discuss the use of 
horizontal interferometers, including atom-chip or guided-atom devices, which potentially allow the 
use of longer coherence times in order to achieve higher sensitivity. 
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I. INTRODUCTION 

Local Lorentz invariance (LLI) in the gravitational in- 
teraction can be viewed as a prediction of the theory of 
general relativity. And it is not a trivial consequence, 
given that alternative theories of gravity have been put 
forward that do not lead to LLI, yet agree with general 
relativity in their predictions for the red-shift, perihelion 
shift, and time delay. Experimental tests of the LLI in 
gravity are required to decide between these theories [l| . 

Another reason to perform tests of LLI in gravity is 
connected to one of the outstanding problems in physics, 
to find a unified theory of quantum gravity. The natural 
energy scale for such a theory is the Planck scale of about 
10^^ GeV. Direct experimentation at the Planck scale is, 
unfortunately, not possible. However, it is possible to 
search for suppressed effects at attainable energy scales 
in experiments of outstanding sensitivity. Violations of 
LLI ("Lorentz violation") are among the relatively few 
candidates for such observable consequences of quantum 
gravity 0. 

The LLI of the non-gravitational standard model has 
been tested for various particles, including photons [1,0, 



electrons [i0Ji,i,[Tl, [nl, [13, protons a[TS [pi, 
neutrons [1, [^, [lOl. Il4l|. and others (see [1, [H, Ha. Il7l.ll8j 

for recent reviews). However, the LLI of gravity itself 
has been studied little to this date. Thus, it remains 
interesting and important to explore the validity of LLI 
in gravity. 

For the theoretical descriptions of the consequences of 
Lorentz violation in gravity in the solar system, a post- 
Newtonian approximation is justified. Such descriptions 
are known as Nordtvedt's anisotropic universe model [Toj 
or th e sta ndard model extension (SME) 0,[2l|,|23], see 
Sec. IIIDI for details. These frameworks use the follow- 
ing Lagrangian to describe the interaction between two 
point-masses M and m (2^ 
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For simplicity, we have taken M to be at rest. We denote 
r the separation between M and m, pointing towards 
TO. The indices j, k denote the spatial coordinates, v 
the relative velocity, and f = r/r. The components of 
s^^ = s'^'^ specify Lorentz violation in gravity. If they 
vanish, LLI is valid. 

The relative weakness of gravity means that only a 
small set of exceptionally sensitive experiments can place 
interesting limits on s. Nordtvedt and Will noted that 



Lorentz violation in gravity would cause a modulation 
of the apparent local gravitational acceleration q as the 
Earth rotates in space. Using gravimeter data [23| taken 
during the international geophysical year, July 1957- 
Dccember 1958, they obtained a limit of \s'^^\ < 4 x 10~^ 
(J, K G {X, Y, Z}) on the spatial components jEdHl- 
More recently, Battat et al. analyzed 34 years of lunar 
laser ranging data, finding bounds on two linear combi- 
nations of the s-"^ at a level of 10^ i° [H]. They also 
found bounds on the three s^'^ at levels of ~ 10~^. 

This means that for four degrees of freedom of s'^^ , 
there has been no improvement for about five decades. 
Moreover, a fundamental issue remains unaddressed: 
testing the isotropy of any force of nature means com- 
paring it to another 'standard' that is assumed to be 
isotropic. In other words, the isotropy of gravity can only 
be tested against the isotropy of another phenomenon, 
such as the velocity of light c. Past tests, however, are 
not easy to analyze in such depth as to make this com- 
parison explicit. For example, it would be very difficult 
to analyze an influence of the isotropy of c on the spring 
gravimeters used by Nordvedt and Will, although a de- 
scription of Lorentz violation in solids is possible in prin- 
ciple [1^, ll^, nil . Hence, it is desirable to perform an ex- 
periment that relies on sufficiently simple physical princi- 
ples so that all channels of infiuence for Lorentz violation 
can be theoretically described. 

With atom interferometry [2^ [s^, HH, we have the 
unique situation of an experiment which is not only sen- 
sitive enough to obtain improved bounds on several ele- 
ments of s'^^, but also 'clean' in the sense that we can 
quantitatively understand the infiuence of Lorentz viola- 
tion in all relevant sectors. The basic principle is shown 
in Fig. [1] A matter wave packet is split by a beam splitter 
to form two interferometer arms. Because of the interac- 
tion with external potentials, such as gravity, the wave 
packet picks up a phase difference Lp between the arms. 
When the arms are recombined at a final beam splitter, Lp 
determines the probability that the atom is found emerg- 
ing from one of the two interferometer outputs. Gravity's 
contribution to can be several 10^ radians for the most 
sensitive devices [H, [s^ . 

Light-pulse atom interferometers use standing waves of 
laser light as gratings to diffract the matter waves; since 
the period of these gratings is given by the laser wave- 
length, the extremely high accuracy of laser frequency 
and phase stabilization can be applied to the measure- 
ment of gravity and other inertial forces. 

Atom interferometers have been used for measure- 
ments of the fine-structure constant a f33.[35l . [3^ . [stI . [ssj , 
the local gravitational acceleration g [33|, its gradient 
[13], the Sagnac effect [MEIil, or Newton's gravita- 
tional constant j4^, i4j] . They rival or exceed the perfor- 
mance of other state-of-the-art methods. Thus, it seems 
natural to apply atom interferometry as a device to probe 
the weakest of the known forces of nature [lOl . Detailed 
proposals have already been put forward, including tests 
of the Einstein Equivalence principle [46[ and detection 
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FIG. 1: Space-time diagram of a Mach-Zehnder atom inter- 
ferometer, t denotes time, z the vertical coordinate. 

of gravitational waves [i^l . 

Here, we will discuss the first atom-interferometer test 
of post-Newtonian gravity [s^. It is also one of the few 
laboratory tests of gravity whose sensitivity is competi- 
tive with the best astrophysics data. This paper is orga- 
nized as follows: We start with a description of the exper- 
iment in Sec. [Hi details not covered there will be found 
elsewhere 33, 39, 48, 49,^. Sections [lEl through |TTC] 
describe the setup; section lTl DI describes the hypothetical 
signal for Lorentz violations in this experiment, including 
the infiuence of the electromagnetic sector. Section III El 
gives details of the data analysis. In Sect. IIIIl we consider 
possibilities for future tests based on horizontal interfer- 
ometers that may be based on atom chips [sil [s^ . [5^ or 
other matter waveguides [s^ [ssj . 

II. EXPERIMENT: A VERTICAL ATOMIC 
FOUNTAIN INTERFEROMETER 

A. Principle 

In our Mach-Zehnder atom interferometer (Fig. [T]), 
two vertical and antiparallel laser beams make the beam 
splitters and mirrors of the interferometer. They drive 
two photon Raman transitions: the atom absorbs a pho- 
ton from one beam and is stimulated to emit a photon 
into the other one. Thus, the atom changes its hyperfine 
state and receives the momentum 2hk of two photons. 
A '7r/2' pulse, which has its intensity and duration cho- 
sen such that this process happens with a probability of 
1/2, acts as a beam splitter; a 'tt' pulse, with a transition 
probability close to one, makes a mirror. 

For our Mach-Zehnder interferometer, a first 7r/2 laser 
pulse transfers the atoms into a superposition of the F — 
3,mF — and F' — A^m'p — states. These move 
vertically relative to each other because of the momentum 
transferred by the laser radiation. A total of three light 
pulses split, refiect, and then recombine the paths to form 
an interferometer. They are separated in time by the 
pulse separation time Tp. 
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The matter waves in both paths acquire a relative 
phase difference if = Sci/h + (see, e.g., Ref. ^] 
for details). The phase of the free evolution of the wave 
packet between the beam splitters is given by the classi- 
cal action Sci- If we restrict our attention to a constant 
gravitational acceleration gi, 



Vibration isolation 



Sci 



2T„ 



1 



mgz 



dt. 



(2) 



z is the vertical coordinate and m the atom's mass. By 
calculating the integral over the classical trajectories, it 
can be shown that this is the same for the upper and lower 
trajectory. Thus, it does not contribute to the phase 
difference between the interferometer arms. 

The phase ipi is because whenever the atom changes 
state during an interaction, the phase of the atom 
changes by an amount equal to the phase of the light 
field. This adds phase when the atom absorbs a photon 
and subtracts phase when the atom emits one. Thus, 
it is easy to see that ipi ^ for the trajectories without 
gravity, depicted by light lines in Fig. [TJ Gravity acceler- 
ates the atoms downwards according to Az — —gP/2, as 
depicted by heavy lines. This gives rise to a phase shift 
of fceffAz, where kes denotes the effective wavenumber. 
To high accuracy, the laser beams can be modeled as 
plane waves, which for two-photon transitions results in 
an effective wavenumber of fccff — ki + ^2, where fci^2 
are the wavenumbers of both lasers. For the upper tra- 
jectory, a photon pair is absorbed at i = and emitted 
a.t t — Tp. We can set the phase of that interaction to 
zero by definition. The phase of the second one will then 
be kesgTp/2. For the lower path, a pair is absorbed at 
t = Tp and one emitted at 2Tp. They add to a phase of 
-k^sgT^/2 + kesg{2Tp)^/2 = 3kcffgT^/2. The difference 
for the two paths is thus (fj — kcsgT^. For the total 
phase shift, we take into account the possibility that the 
phases of the laser pulses at z = can experimentally be 
set arbitrary values tpi, iP2 , by suitable control of the 
laser system. We obtain [S^ 



= kcsgTp - 

= ipi- 2ip2 



+ <P3- 



(3) 



For our experiment, Tp = 0.4s and k ~ 27r/(852nm), so 
~ 2.3x 10''rad. In the experiment, we set ipL = rT^ by 
ramping the difference frequencies of the laser (actually, 
the continuous ramp can be approximated by discrete 
steps, as the laser pulse duration is very short). To mea- 
sure g, the interferometer phase ip is zeroed by adjusting 
the ramp rate to ro « 27r x 23MHz/s. This corresponds 
to finding the center of the interference pattern. Then, 
g — ro/(fccff). If ip can be measured to, e.g., 1 mrad, we 
obtain a resolution of 10"^" in g. 



B. Fountain 

In our experimental setup (Fig. [5]), we assemble about 
10^ Cs atoms within 650 ms from a background vapor 
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FIG. 2: Setup. A/4; 1/4- wave retardation plate. 
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pressure of ^ 10^^ mbar in a 3-dimensional magneto- 
optical trap (3D-M0T). A moving optical molasses 
launch accelerates them vertically upwards to a 1— s 
ballistic trajectory with a temperature of 1.2-2 /iK. Ra- 
man sideband cooling in a co-moving optical lattice re- 
sults in ~ 3 X 10^ atoms in the F = 3, mp — 3 state at a 
(3D) temperature of 150 nK that form a cloud of roughly 
3mm^ area \0l. A sudden change in the magnetic field 
followed by a 120-/is microwave pulse transfers ^ 20% of 
them into the F = 4, tof = state. Atoms left over in 
the F = i state are then cleared away using a resonant 
laser pulse. A solenoid generates a small magnetic bias 
field to set the quantization axis. 



C. Interferometer 

The off-resonant Raman pulses for the beam splitters 
have a wavelength of 852 nm and are generated by two 
grating-stabilized extended cavity diode lasers that are 
based on 100-mW laser diodes SDL-5411. The first is 
frequency stabilized ('locked') to a cesium vapor cell us- 
ing Doppler-free saturation spectroscopy. It arrives at 
the experiment with a detuning of -1030MIIz from the 
&Si/2, F 6P3/2, F' = A transition in Cs. The sec- 

ond one is phase locked to the first one with a frequency 
difference close to the hyperfine splitting of ~ 9192 MHz, 
referenced to a LORAN-C frequency standard. The two 
lasers are overlapped on a beam splitter and the com- 
bined beam, containing 20 mW of each one, is trans- 
mitted to the experiment via a single-mode, polarization 
maintaining optical fiber. There, the beams are switched 
and intensity-controlled by an acousto-optical modulator 
(Isomet 1205). They are coUimated with a 1/e^ intensity 
diameter of about 2.5 cm, and pass the vacuum cham- 
ber with linear polarization. Retro-reflection on top of it 
with two passes through a quarter- wave retardation plate 
forms a lin_Llin polarized counterpropagation geometry 

(Fig. m- 

The error signal for the phase lock is generated by de- 
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FIG. 3: Typical fringe obtained in our experiment 

tecting the beat note after the fiber. Phase noise sources 
before this point can thus be taken out by the phase lock 
loop, whereas most others are common to both beams 
and do not affect the interferometer. The most impor- 
tant exception are residual vibrations of the top mirror, 
which we therefore reduce to below 5 x 10^^(?/\/Hz in 
a frequency range of 0.1-lOHz by a sophisticated active 
vibration isolator [isj . 

For fluorescence detection with a Hamamatsu R943-02 
photomultiplier tube (PMT), the F = 4 interferometer 
output is driven on a cycling transition; the F = 3 out- 
put is detected subsequently after optical pumping to the 
F = 4 state. Normalization of the signals takes out vari- 
ations in the number of launched atoms. 

Fig. [3] shows a typical gravity fringe with a pulse 
separation time of Tp = 400 ms, taken with 40 launches 
that take 75 s total. The sinewave-fit has a phase uncer- 
tainty of 0.031 rad, and determines g to an uncertainty 
of ^ 1.3 X 10~^g. This corresponds to 11 x lO^^g/VHz. 
An improved short-term resolution of 8 x lO^^g/v^Hz 
can be reached by taking data at the 50% points of the 
fringes only. However, as this method is more sensitive 
to systematic effects such as drift of the PMT sensitivity 
[39| . we used fringe- fitting for taking the long term 
data. Our resolution is about three times better than 
the best previous one, which was reported by Peters et 
al. [3§|. This is mainly a consequence of the increased 
interaction time T — 400 ms. It also surpasses the best 
classical absolute gravimeter, the FG-5 falling corner 
cube gravimeter, by a factor of about 20. 



D. Hypothetical Signal 

The notion that the gravitational force between two 
objects might depend on the direction of the separation 
r could be described in very simple terms. For this work, 
however, we want to use a model that is as general as 
possible on the one hand and compatible with accepted 
principles that underlie the standard model and gravita- 



tional theory on the other hand. Two such models sug- 
gest themselves, Nordtvedt's anisotropic universe model 
[iil and the standard model extension (SME) [23,[lll[2l| 
(a simplified formalism with one parameter has also been 
proposed [1^.) The SME starts from a Lagrangian for- 
mulation of the standard model and gravity, adding gen- 
eral Lorentz violating terms that can be formed from the 
fields and tensors. It is therefore a comprehensive the- 
oretical framework that allows us to model violations of 
LLI in the various sectors of the standard model and 
gravity 

As Lp — kcsgTp, the hypothetical signal for violations 
of LLI can be due to a change in g and fcoff. The sig- 
nal may thus be due to Lorentz violation in both the 
electromagnetic as well as the gravitational sector. 



1. Gravitational sector 

In a post-Newtonian approximation, the Lagrangian 
for the gravitational interaction in the SME is given by 
Eq. [TJ The two-body Lagrangian of the anisotropic uni- 
verse model is similar, but = and the coefficients of 
vl and ^abW^f^j are independent of each other. 

In principle, the components of s can be defined in any 
inertial frame of reference. For experiments on Earth (as 
well as on satellites), it is convenient to choose a Sun- 
centered celestial equatorial reference frame A]. It has 
the X axis pointing towards the vernal equinox (spring 
point) at Oh right ascension and 0° declination, the Z 
axis pointing towards the celestial north pole (90° dec- 
lination) and the Y axis in the way needed to complete 
the right handed orthogonal dreibein. Earth's equatorial 
plane is the X — Y plane and the orbital plane of the 
Earth is tilted at an angle i] ~ 23° with respect to the 
latter. The time scale T is set by T = when the Sun 
passes the spring point, which, for example, happened 
in 2001 on March 20, 13 h 31 min universal time (UT). 
Sun-centered frame quantities have Greek (0-3) or capi- 
tal Latin indices. 

We also define a laboratory frame, which has the x axis 
pointing south, the y axis east, and the z axis vertically 
upwards. The laboratory time scale is set by = at 
any one instant when the y and the Y axis coincide. The 
difference between the two time scales and T can 
be written as a phase difference = ojq{Tq — T) [2^ . 
where 0^0 is the sidereal angular frequency of Earth's 
rotation; (f> ~ —1.77 for this experiment. Laboratory- 
frame quantities have small Latin indices. 

The derivation of the time-dependent modulations of 
g for an observer on Earth involves taking into account 
the rotation and orbit of the Earth; the Earth itself is 
modelled as a massive sphere having a spherical moment 
of inertia of J® « [13] (not to be confused with 

the conventional moment of inertia, which for Earth is 
about Mi^r'^/i). It suffices to consider the first order 
in the Earth's orbital velocity ~ 10~*. Bailey and 
Kostelecky ;22] have studied this in detail, and we refer 
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the reader to this reference for the detailed signal com- 
ponents in the purely gravitational sector. 



2. Electromagnetic sector 

To study the variations of fc^ff caused by Lorentz vio- 
lation, we start from the Lagrangian for the electromag- 
netic sector of the SME, 



C 



1 



(4) 



where F^" is the electromagnetic field tensor and the 
vector potential. The second term is proportional to a 
dimensionless tensor {kp)K\f^i,, which vanishes, if Lorentz 
invariance holds on electrodynamics. The tensor has 19 
independent components A brief summary of the 
Maxwell equations that are derived from this Lagrangian 
is given in appendix [XI 

For studying the Lorentz-violating modification to the 
effective wavevector k^ff in the atom interferometer, the 
plane wave solutions have to be found. Making the stan- 
dard ansatz F^v{x) — i^^,y(p)e^*'^°^ for a plane wave 
with a wave 4- vector fc" — (fc*^, k) and inserting into Eq. 
(|A1[) one obtains the dispersion relation, that will allow 
us to determine k^s as a function of the direction of prop- 
agation. Let 



= {kFr^f^-p^p,, r^^- 

m 

Then the dispersion relation is Q 
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il + p±a)\k\ 



(5) 



(6) 



The last term in this relation, which is proportional to 
a, is purely polarization-dependent. Astrophysics shows 
that such a dependence, if it exists, is well below the 
levels relevant here |3:j. We can thus assume a = 0. 

To obtain the explicit time-dependence of the wavevec- 
tor in our experiment, we need to transform the quanti- 
ties, which are conventionally defined in a sun-centered 
celestial equatorial reference frame frame into the labo- 
ratory frame [1]. 

Explicit values of the Lorentz violating quantities de- 
pend on the definition of coordinates and fields [1, 0, 
El, ■ The freedom to define these can be used to set 
certain components to zero. A particular definition, for 
example could be made by requiring K^i^ = (see ap- 
pendix[X| . A choice like this would in general make these 
terms reappear in other sectors, such as the gravitational 
and fermionic sectors. In the following, we do not make 
any particular assumptions on such definitions and retain 



TABLE L Signal components for vertical atom interferome- 
ters. X is geographical colatitude, 42.3° for Stanford. 
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3. Combined signal 

Adding the contributions of the electromagnetic and 
the gravitational sector, the time-dependence of the in- 
terferometer phase can be expressed as a Fourier series 
for the time-dependence [23] 



^ C,„ cos(w„jt ^- 0,„) -f D„i sin(u;„j< -I- 0,„). (7) 



The coefficients C„i , for the six frequencies m S 
2^0, ± rj, ± rj} are functions of the Lorentz 
violations, that are given by the components of s^^ and 
{kp^KXyiv and the frequencies of Earth's orbit fi^ = 
27r/(ly) and rotation ~ 27r/ (23.93 h). For a vertical 
interferometer, we obtain the signals for Lorentz viola- 
tion in electromagnetism and gravity, see Tab. [H We 
denoted 14 = 1 - 3/®/(M0i?|) « -1/2. The symbols 
Ke_ , Ko_|_ denote linear combinations of the elements of 
{kF)K\tiv that are defined in appendix [XI 

It turns out that for these components (but not, for 
example, for Co, ZJq or the ones in Tab. IVip . the substi- 
tutions 
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all the quantities s and Kg-, 



in full generality. 



(8) 

can be used to obtain the combined signal components 
from the purely gravitational ones listed in Tab. IV of 
Ref. (H. 



E. Data analysis 

The combined data spans about 1500 d, but frag- 
mented into three relatively short segments: ^60 h of 
data taken with this setup, as well as a ~ 60h and a 
'--^lOd run reported previously [1^, see Fig. U) 

Major systematic effects in this experiment are tidal 
variations of the local gravitational acceleration. Sub- 
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FIG. 4: Data in 10 ^g. Each point represents a 60-s scan of one fringe (75-s after TJD1900). 



traction of a Newtonian model of these influences [58| 
based on the relative positions of the Sun, the Moon, 
and the planets, yields the residues shown at the bot- 
tom of Fig. m In addition to the tidal model used in 
our Letter, we here use an additional model of the local 
tides [59|. This will allow us to significantly reduce the 
relatively large estimates on some of the a coefficients 
reported in the Letter. As our old analysis and analy- 
ses using lunar laser ranging [25l| , the present analysis is 
based on an assumption of no accidental cancelations of 
signals, i.e., that the model does not contain the same 
influence of Lorentz violation as our measurement. 

In order to extract the signal components of Tab. HI 
we calculate the discrete Fourier transform. As is typical 
of a finite set of data, the Fourier components overlap. 
This overlap can be quantified by a covariance matrix. 
For a compact presentation, let 

Um = (2Lj©,a;©, 2we -f ri, 2w0 - ri, -I- fi, W0 - fi) 
a2m+i = cos{ujmt + (t)m.), = sin(ti;„i< + 0„) . (9) 

For a total of K data, let d{tk),k G {l,...K} be the 
datum taken at the time tk- The sine and cosine Fourier 
transforms are 



TABLE II: Correlation matrix cov. 



Z ^ T 



(10) 



the dm that are obtained from the fragmented data are 
a linear combination 
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(11) 



^cov{am,a„)D„. 

n 

of the corresponding Fourier components for hypo- 
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0.6 


0.6 


0. 


0. 


0. 


0. 


0.6 


0.3 


0. 


1 



thetical continuous data. The covariance matrix 



cov(am,a„) = — am{tk)an{tk)j 



(12) 



where the sum is over all tk- For data that span a large 
time-scale, it approaches a unit matrix, i.e., the overlap 
of the Fourier coefficients becomes negligible. The co- 
variance matrix for our data is shown in Tab. |TT1 

Multiplication of the covariance matrix (Tab. |TT]) 
with the signals listed in Tab. [1] (written as a vector 
{C2u], D2lj, Cu, . . . , Duj-q)) gives the linear combinations 
that the experiment is sensitive to, see Tab. IIIIl For 
obtaining bounds on Lorentz violation, we perform a nu- 
merical Fourier analysis of the data for the 12 compo- 
nents Cm-, dm, see the table. The error is estimated by 
performing a Fourier analysis at several frequencies above 
and below the signal frequencies and computing the root 
of the mean square. 

The limits listed in Tab. IIIIl are on linear combina- 
tions of parameters. The 12 results are sufficient to de- 
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TABLE III: Signal components as obtained by a numerical Fourier transform of the data. They correspond to linear combina- 
tions of the components of a, with linear coefficients as tabulated. The result of the Fourier transform is listed as fraction of 

a- 



Comp. 


Tx — 


cr-' ^ 




Xa yy 
a ~ a 








Meas. /10~^ 




-0.091 


0.040 


- 0.002 


- 0.079 


0.004 


- 0.009 


0.004 


-0.097 ±0.24 




-0.039 


- 0.094 


- 0.007 


0.001 


- 0.241 


- 0.007 


0.001 


0.009 ± 0.24 




-0.066 


0.001 


- 0.142 


- 0.003 


- 0.007 


- 0.240 


0.001 


0.69 ±0.27 




0.003 


- 0.063 


- 0.136 


0.001 


0.001 


0.001 


- 0.243 


-0.034 ±0.27 




-0.105 


- 0.094 


- 0.005 


- 0.020 


- 0.141 


- 0.007 


- 0.004 


0.19 ±0.24 


d2uj+n 


0.098 


- 0.107 


0.011 


0.044 


- 0.068 


0.016 


- 0.003 


-0.22 ± 0.24 


C2u}-n 


-0.005 


0.158 


- 0.001 


- 0.021 


0.144 


- 0.002 


0.006 


-0.15 ±0.24 


d2uj-n 


-0.157 


- 0.000 


- 0.004 


- 0.043 


- 0.064 


- 0.014 


0.006 


0.40 ± 0.24 




-0.065 


- 0.034 


- 0.131 


- 0.003 


- 0.002 


- 0.140 


- 0.143 


0.20 ±0.27 


d^+n 


0.036 


- 0.064 


- 0.053 


0.003 


0.004 


0.139 


- 0.151 


-0.73 ±0.27 


Cw-n 


-0.067 


0.029 


- 0.003 


- 0.005 


- 0.001 


- 0.143 


0.140 


0.28 ±0.27 




-0.029 


- 0.061 


- 0.231 


- 0.001 


- 0.007 


- 0.142 


- 0.148 


0.25 ±0.27 



termine all 7 parameters. In order to obtain indepen- 
dent estimates for the parameters while making opti- 
mum use of the experimental data, we proceed as fol- 
lows. Each datum corresponds to a probability dis- 
tribution pnia"^^ , . . . , cr^^), where n = 1, ... 12, that 
we assume to be Gaussian with the center and stan- 
dard deviation as tabulated. Multiplying all 12 dis- 
tributions results in one overall probability distribution 
P{(T^^ , . . . ,a^^) = AYll^^iPn, where A is a normal- 
ization factor. In order, for example, to obtain an es- 
timate on a'^-^ that is independent of the other coeffi- 
cients, we integrate over the other variables, P{a^^) — 



J — oo J —oc ^ 



a ) (where all 



variables but a^-^ are integrated over and A' is another 
normalization factor). The result is a probability dis- 
tribution for a^^ that has one maximum and standard 
deviation, which are our estimate and error bar for this 
coefficient. Specifically, since we assumed the p„ to be 
Gaussian, P{a^^) will be a Gaussian and the most prob- 
able estimate as well as one a error are read off. 

Tab. IIVI shows the results thus obtained. Compared 
to the ones that we published in our Letter [SS^ (derived 
from the same data), the limits are now more uniform: 
whereas the best limits on the three a^"^ and four (com- 
binations of) o"^^ are of the same order, the worst ones 
from the new analysis are about four times better than in 
the Letter. This is a result of our use of a more sophisti- 
cated tidal model as well as an optimum method to derive 
individual limits from the limits on linear combinations. 



F. Combination with lunar laser ranging (LLR) 

Our experiment thus provides the only measurement 
of the components of a so far; these are combinations 
of gravitational and electromagnetic Lorentz violation. 
To compare our data with previous experiments, which 
have not been analyzed for the electromagnetic influence 
(although such an influence exist), we here adopt the 
assumption that there is no Lorentz violation in electro- 



TABLE IV: 



Coeff. 


^TX 


(-3.1 ±5.1) X 10"^ 




(0.1 ±5.4) X 10"^ 


^TZ 


(1.4 ±6.6) X 10"^ 




(4.4 ± 11) X 10"^ 




(0.2 ±3.9) X 10"^ 




(-2.6 ±4.4) X 10"^ 




(-0.3 ±4.5) X lO"'' 



magnetism, i.e., k'l_ — k'^^ = 0. Our limits (Tab. IIV| 
then correspond directl y to bounds on s. 

Lunar laser ranging |25j bounds two linear combina- 
tions of s-'^ (Hi 



^22 



-XX 



O.U8(s"" + s 

-XX 



YY 



-0.31(s---s-^^)-1.7s^^ 
+0.60s-^^ + 0.42s^^ 
(1.3 ±0.9) X 10^1°, 
0.43(s^^ - s^^) - 0.31s^ 
-0.23s^^ - 0.33s^^ 
(6.9 ±4.5) X 10"". 



7;XY 



(13) 



They can be combined with our results from Tab. IIIII 
(using the method described above) to obtain indepen- 
dent limits on one more degree of freedom of s"'^ . Battat 
et al. also report four limits on the three s'^'^, [25| 

^ -0.60s^^ ± 0.82s'^^ 

= (-0.8 ± 1.1) X 10-6, 

= -0.535^^-0.755^^ 

= (-5.2 ±4.8) X 10-^ 

sn^c = -3.1s^^- l.Ls^^ 

= (0.2 ±3.9) X 10"^ 

SQ^s = -3.4s^^ ± 0.037s^^ 

= (-1.3 ±4.1) X lO-'^. (14) 



0.40s 



TZ 



0.094s 



TX 



0.15s 



TZ 
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TABLE V: Bounds resulting from combining our data with 
the ones from lunar laser ranging as reported by Battat et al. 
[25j . assuming vanishing Lorentz violation in electrodynamics. 



Coeff. 




(0.5 ±6.2) X 10" '' 




(0.1 ± 1.3) X 10"'' 


-TZ 


(-0.4 ±3.8) X 10"'' 




(-1.2 ± 1.6) X 10"^ 


-XX ^ -YY _ 


(1.8 ±38) X 10"^ 




(-0.6 ± 1.5) X 10"^ 




(-2.7 ± 1.4) X 10"^ 




(0.6 ± 1.4) X lO"'' 



They can be combined with ours to increase the resolu- 
tion of the hmits. 

Tab. |V] lists the results thus obtained. They repre- 
sent the most complete bounds on Lorentz violation in 
gravity, providing individual limits on the s as well as 
more components of s and higher resolution than either 
experiment. The only degrees of freedom of s''^ that are 
not bounded are s'^'^ and the trace, which do not lead to 
signals to first order in the Earth's orbital velocity. 



III. SIGNAL FOR HORIZONTAL 
INTERFEROMETERS 

In this section, we consider tests of gravity with hor- 
izontal atom interferometers, including guided atom de- 
vices. Testing LLI in gravity is a task that makes good 
use of the features of such interferometers, in particular 
long coherence times and hence high resolution. More- 
over, since the signal for violations is a time-dependent 
modulation, the stability of the interferometer on time 
scales much larger than the modulation frequencies are 
not a primary concern for such tests. 

For simplicity, we shall again assume a vanishing of 
Lorentz violation in electrodynamics throughout this sec- 
tion. 

A test of the LLI of gravity can be performed by mea- 
suring a Lorentz-violating horizontal acceleration. These 
accelerations are given by [2^ 1 



R(B sin X cos X 



■Tz 



■ -Tx-irz 



(15) 



where 13 = 1 — /^/(M^i?^) w 1/2. For the purpose of 
this section, we can take = g as well as a = \a\ Kg 
to be constant. Such tests can, for example, be based 
on a torsion pendulum, which is suspended off its center 
of mass. Nevertheless, they might reach superior sensi- 
tivity compared to vertical gravimeters such as the atom 
interferometer discussed previously. 

Compared to a conventional torsion pendulum, such 
experiments involve special challenges associated with 



maintaining the pendulum within the horizontal plane. 
However, measurement of the horizontal accelerations 
that are given by Eqs. (jlSp with atom interferometry 
is possible using conventional horizontal interferometers. 
Moreover, interferometers in the horizontal plane can 
be built well using atom-chip or atomic waveguide tech- 
niques. In contrast to atomic fountains, they allow long 
pulse separation times Tp in a compact setup. It is there- 
fore interesting to study the signals for Lorentz violation 
in post-Newtonian gravity and electromagnetism for such 
an interferometer. 

We assume a horizontal Mach-Zehnder interferometer 
with the laser beams pointing into a direction of 



(cos 9, sin 6*, 0) 



(16) 



in the laboratory frame. As before, the phase shift is 
given by k^sTp {x ■ g), where the local gravitational accel- 
eration g has vertical as well as horizontal components. 

The calculation of the induced time-dependence of the 
interferometer phase proceeds via the transformations 
between the laboratory frame and the sun-centered stan- 
dard frame. The fastest way to do this is probably by 
analogy to the case of torsion balances that has been 
considered in [24I; see appendix IB] After all, both mea- 
sure the accelerations given by Eqs. As a result, 
we can express the contribution of Lorentz violation in 
gravity to the phase as 

ifi = kcshgTp ^ {[En sin a„ - En cos a„] sin a;„T 

n 

— [En cosa„ ± En sina„] cosa;„T) . (17) 

The amplitudes and phases in this expression are given 
in Tab. IVII It is evident that horizontal interferometers 
provide access to four independent linear combinations 



of 



-JK 



(the same ones as vertical interferometers) and 



sufhcient data to determine all the s^"^. 



IV. SUMMARY AND OUTLOOK 

In this paper, we have presented atom-interferometry 
tests of the local Lorentz invariance of post-Newtonian 
gravity and electrodynamics. As a comprehensive and 
quantitative model for violations, we use the standard 
model extension [20, 21]. The relevant violations of LLI 
in gravity are encoded by a tensor s (2^ ; those for viola- 
tions in electromagnetism are expressed by the matrices 
Ke_ and Ko+ 3. The experimental signal for Lorentz 
violation is a time-dependence of the local gravitational 
acceleration as the Earth orbits in the solar system. We 
discuss an experiment that has been performed by us (ssj 
as well as possible future experiments with horizontal in- 
terferometer geometries. 

Our experiment is a vertical Mach-Zehnder atom inter- 
ferometer. It uses a bright source of cesium atoms in a 
1-m high atomic fountain and a pulse separation time of 
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TABLE VI; Amplitude and phase for the gravitational part of the signal in horizontal interferometers 



Amplitude 

'eI^ -s'^'' sine sin X- 3(5"^'^ - s*' ) cos 61 sin 2x 2^ 

F2^ \ (s-^-^ - -s^^) sin e sin X - ^s^^ cos Q sin 2x 

Eu — s^'^ sin ^cos X — s'^'^ cos cos 2x (p 

Fuj s"^^ sin S cos X — s'^^ cos cos 2x 4> 



E2Lo+n ^V^s^"^ sin6l(l — COS77) sinx — ^V'es^^ cos6l(l — cost;) sin2x 20 

F2ui+n I V®s^^ sin 61(1 — cos 77) sin x + 3 VeS^"^ cos 61(1 — cos rj) sin 2x 20 

E2uj-n — |V®s^^ sin6l(l + cos»7) sinx + 3:V©s^-*' cos6'(l + COST)) sin2x 20 

F2ui-n — iV®s^^ sinS(l + cosr/) sinx — ^Ves^"^ cos6'(l + COS77) sin2x 20 

Euj+n — 5 Ves^^ sin 6^ sin r) cos x + Ves^^ cos 6^ sin 77(i — cos^ x) 4> 
— i V^s^'^ sin^(cos 77 — 1) cos x 

-FL+n i Ves^"^ sin&sin?7cosx — V®s^^ cos ^(1 — cos 77)(| — cos^ x) <^ 

+V(5S^^ cos sin 77(1 — cos^ x) 

E^-n — iV®s^^ sin 6(1 + COS?;) cosx + ^Ves^^ sinSsin J7C0SX 

+Ves^'^ cosSsin77(i — cos^ x) 

Fui-n iVes^'^ sin6'sin77cosx + Ves"^ cosSsin77(i — cos^x) 
+V©s'^^ cose(l + cosj7)(i -cos^ X) 

En — \ V® s"^^ cos 9 cos 77 sin 2x + V® s ^ cos 6* sin 77 sin 2x 

Fn ^t/es^-^ cose sin 2x 0_ 



400 ms. A resolution of up to 8 x 10^^ g/VHz is reached, 
the highest of any cold-atom gravimeter so far. 

For this experiment, the signal for Lorentz violation is 
given by a particular combination of coefficients entering 
the gravitational and electromagnetic sectors: o"^^ = 
s-'^ - ji^ and a^-' = s^-^ + ejKLKo+ /{2i4). Here, 
14 ~ —1/2 is given by the Earth's spherical moment of 
inertia. 

For the analysis, we use about 2.5 d of data taken with 
the setup just described as well as about 12.5 d of data 
taken with a previous setup (39j . The data have been 
taken over a total time interval of about 1500 d. A ma- 
jor systematic effect are the tides, which we account for 
by subtracting a Newtonian model that is based on the 
relative positions of the Earth, the Moon, the Sun, and 
planets. By Fourier analyzing the residuals at different 
combinations of the frequencies of Earth's rotation and 
orbit, we obtain bounds on twelve linear combinations of 
seven components of o"^^ at the 10~^ level. 

Our limits are the strongest bounds on several combi- 
nations of coefficients for Lorentz violation, even in view 
of p revious work from geophysical gravity observations 
[l9| and lunar laser ranging [25|. This makes our experi- 
ment one of the few competitive laboratory tests of grav- 
ity. Moreover, our experiment is presently the only test 
of Lorentz invariance of gravity that explicitly takes into 
account the possibility of Lorentz violation in the non- 
gravitational standard model. This is important when 
interpreting the results in the context of the search for 
quantum-gravity signals. 

The physics behind atom interferometry and lunar 
laser ranging is similar, monitoring the trajectory of a 
proof mass within Earth's gravitational field by laser 
beams. However, they differ in the orbit (if one can 
think of the atoms' trajectory as an orbit) and quantum- 
mechanical nature of the mass. Indeed, the equivalence 



principle for quantum objects has been discussed [60|. 
Also, ours is a laboratory experiment, which typically 
offers superior control over the experimental systemat- 
ics. Moreover, it is so far the only experiment where the 
simultaneous influence of the non-gravitational and grav- 
itational effects are understood quantitatively and which 
accordingly states combined bounds. 

Assuming that Lorentz invariance in electrodynamics 
vanishes, the data from LLR and our experiment can be 
combined to yield improved and more detailed bounds. 

We also studied the signals for Lorentz violation (in 
gravity and electromagnetism) for horizontal interferom- 
eters. They are attractive for this type of measurement, 
as they offer increased coherence time in a compact setup. 

Taking several months worth of data would help to 
eliminate the dominant lunar tides, the major systematic 
effect. With a dataset that spans a year, even the leading 
solar tides could be suppressed. 

Future work may lead to alternative bounds from tor- 
sion balances, gravimetry data that is routinely taken in 
geophysical research, or an analysis of data from the grav- 
ity probe-B satellite [6l|. Note, however, that our results 
are still not limited by the any fundamental influence 
such as quantum projection noise. With typically 10^ 
atoms per launch, a quantum projection limited gravime- 
ter could reach the 10~^^g level per launch and 10~^'^g 
per day, if other noise sources (notably phase noise and 
vibrations) can be controlled. Also, the performance of 
atom interferometers can be increased further by utiliz- 
ing beam splitters that transfer the momentum of many 
photons, thus leading to sensitivity increases by factors of 
12 (demonstrated) to 100 times (anticipated) [H, HI, HI] . 
This promises improved tests of gravity based on atom 
interferometry, deepening our understanding of the fun- 
damental principles of Nature. 
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APPENDIX A: ELECTROMAGNETIC SECTOR 

The Maxwell equations in vacuum that are derived 
from the Lagrangian, Eq. are 

d^F'^'^ - 0, (Al) 

where 

pf"" = ie'^^^'^F^^. (A2) 

They can be written in analogy to the Maxwell equations 
in anisotropic media ^] : With the 3x3 matrices 

(latin indices take the values 1,2,3) one can define D and 
H fields 

fD\_fl + KDE KDB \fE\ 

[h )^[ ^HE t + ^^HB )[b)' ^^'^> 

where 1 represents the 3x3 unit matrix. The Maxwell 
equations can now be expressed as Q 

V X H - doD ^ 0, V • B = 0, 

V X E + doB ^0, V • = 0. 

Lorentz violation in electrodynamics is thus analogous to 
electrodynamics in anisotropic media. For later use, we 
define the linear combinations 





1 , 

-^\*^DE + 


KHBy'', 




1 , 

-^\*^DB + 


K.HEy*' , 




1. 

-^\*^DE - 




if^o-y = 


\{>^DB - 


HHEy'', 


Ktr — 


\{^DEf 





Of these, the ten degrees of freedom of iio- and Ke+ en- 
code birefringence; they are bounded to below 10"'^'' by 
observations of gamma-ray bursts [3, The residual 
nine cause a dependence of the velocity of light on the 
direction of propagation or, vice versa, such a depen- 
dence of the wavenumber of light having a constant fre- 
quency. They are therefore relevant in interferometry 
experiments. 

APPENDIX B: ANALOGY OF HORIZONTAL 
INTERFEROMETERS AND TORSION 
PENDULUM EXPERIMENTS 

The free oscillations of a torsion pendulum are gov- 
erned by the differential equation "i? -I- {k/ I)d — t, where 
/ is the moment of inertia, and r represents the driving 
force. It can be expressed by an effective potential V as 
r = —dV/M. In the limit of k ^ oo, we simply have 

Kostelecky and Bailey [l^ give the oscillations 'd{T) by 
their Eq. (134). For small excitations i?, 

^ = -mroia'^SN - a^CN + ^acul sin^ xSn2)- (B2) 

Equating this with Bailey and Kostelecky's Eq. (134) in 
the limit of — ^ oo, we obtain 

= *35 X! ( ~ ^" sin(w„T) 

n 

-[En cosa„ -I- Fn sina„] cos(w„T)) (B3) 

where the amplitudes En, Fn and the phases an are tab- 
ulated in Tab. V of [22] ■ They represent the time- 
dependence of the Lorentz-violating horizontal acceler- 
ations as functions of the Earth's orbital parameters and 
the Cm,Sn,S2n- Specializing to the case of one mass. 
Cat = cos 9, Sn = sin 9, and S2N = sin 29. By com- 
parison of coefficients, we obtain explicit expressions for 
ax,a,y, and in the laboratory frame, that we can in- 
sert into the equation ip — kcsgTp for the interferome- 
ter's phase. Note that with our specification of the atom 
interferometer's orientation x, Eq. p6|) . 9 = {} + tt/2, 
which has been taken into account in Tab. IVII 
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